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1 Introduction 

It is twenty years ago almost to the day that the les Houches school on Fields, Strings and 
Critical Phenomena took place. It came on the heels of a frenzied period of five years or 
so following the seminal paper of Belavin, Polyakov and Zamolodchikov (BPZ) in which 
the foundations of conformal field theory and related topics had been laid down, and 
featured lectures on CFT by Paul Ginsparg, Ian Affleck, Jean-Bernard Zuber and myself, 
related courses and talks by Hubert Saleur, Robert Dijkgraaf, Bertrand Duplantier, Sasha 
Polyakov, Daniel Friedan, as well as lectures on several other topics. The list of young 
participants, many of whom have since gone on to make their own important contributions, 
is equally impressive. 

Twenty years later, CFT is an essential item in the toolbox of many theoretical condensed 
matter physicists and string theorists. It has also had a marked impact in mathematics, in 
algebra, geometry and, more recently, probability theory. It is the purpose of these lectures 
to introduce some of these tools. In some ways they are an updated version of those I 
gave in 1988. However, there are some important topics there which, in order to include 
new material, I will omit here, and I would encourage the diligent student to read both 
versions. I should stress that the lectures will largely be about Conformal Field Theory, 
rather than Conformal Field Theories, in the sense that I'll be describing rather generic 
properties of CFT rather than discussing particular examples. That said, I don't want the 
discussion to be too abstract, and in fact I will have a very specific idea about what are the 
kind of CFTs I will be discussing: the scaling limit (in a sense to be described) of critical 
lattice models (either classical, in two dimensions, or quantum in 1+1 dimensions.) This 
will allow us to have, I hope, more concrete notions of the mathematical objects which 
are being discussed in CFT. However there will be no attempt at mathematical rigour. 
Despite the fact that CFT can be developed axiomatically, I think that for this audience 
it is more important to understand the physical origin of the basic ideas. 



2 Scale invar iance and conformal invar iance in criti- 
cal behaviour 

2.1 Scale invar iance 

The prototype lattice model we have at the back of our minds is the ferromagnetic Ising 
model. We take some finite domain T> of d- dimensional euclidean space and impose a 
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regular lattice, say hypercubic, with lattice constant a. With each node r of the lattice 
we associate a binary-valued spin s(r) = ±1. Each configuration {s} carries a relative 
weight W({s}) oc exp(X)rr'ex> ^( r — r')s(r)s(r)) where J{r — r') > is some short-ranged 
interaction (i.e. it vanishes if \r — r'\ is larger than some fixed multiple of a.) 

The usual kinds of local observable 0' at (r) considered in this lattice model are sums of 
products of nearby spins over some region of size 0(a) , for example the local magnetisation 
s(r) itself, or the energy density Xv J{ r ~ r')s(r)s(r'). However it will become clear later 
on that there are other observables, also labelled by a single point r, which are functions 
of the whole configuration {s}. For example, we shall show that this Ising model can 
be mapped onto a gas of non-intersecting loops. An observable then might depend on 
whether a particular loop passes through the given point r. From the point of view of 
CFT, these observables are equally valid objects. 

Correlation functions, that is expectation values of products of local lattice observables, 
are then given by 

<03r*(ri)#*(r 2 ) • • • ^\r n )) = Z^ £ $* (n) . . . ^ (r n )W ({s}) , 

M 

where Z = Y^{s] W({s}) is the partition function. In general, the connected pieces of these 
correlations functions fall off over the same distance scale as the interaction J(r — r'), but, 
close to a critical point, this correlation length £ can become large, ^> a. 

The scaling limit is obtained by taking a — > while keeping £ and the domain T> fixed. 
In general the correlation functions as defined above do not possess a finite scaling limit. 
However, the theory of renormalisation (based on studies in exactly solved models, as well 
as perturbative analysis of cut-off quantum field theory) suggests that in general there 
are particular linear combinations of local lattice observables which are multiplicatively 
renormalisable. That is, the limit 



lim a' E;=i *' ^T( ri )4> l f(r 2 ) . . . ^ (r n )> (1) 

a — »U 

exists for certain values of the {xj}. We usually denote this by 

(0i(ri)0 2 (r 2 )...0 n (r n )>, (2) 

and we often think of it as the expectation value of the product the random variables 
4>j(rj), known as scaling fields (sometimes scaling operators, to be even more confusing) 
with respect to some 'path integral' measure. However it should be stressed that this 
is only an occasionally useful fiction, which ignores all the wonderful subtleties of renor- 
malised field theory. The basic objects of QFT are the correlation functions. The numbers 
{xj} in (CQ) are called the scaling dimensions. 

One important reason why this is not true in general is that the limit in (TjQ) in fact only 
exists if the points {r,} are non-coincident. The correlation functions in (J2J) are singular 
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in the limits when r\ — > Tj. However, the nature of these singularities is prescribed by the 
operator product expansion (OPE) 

(Mr,)U r i) •••> = £ Cijkirt - r 3 )(M(r* + r d )/2) . . .) . (3) 

k 

The main point is that, in the limit when [r*i — 7~j j is much less than the separation between 
Ti and all the other arguments in . . ., the coefficients C^ are independent of what is in 
the dots. For this reason, (J3J) is often written as 

M r i) • <t>j( r i) = J2 C ijk(ri ~ rj)(f) k (( r i + r j)/ 2 ) , ( 4 ) 

k 

although it should be stressed that this is merely a short-hand for (J3j). 

So far we have been talking about how to get a continuum (euclidean) field theory as 
the scaling limit of a lattice model. In general this will be a massive QFT, with a mass 
scale given by the inverse correlation length £ _1 . In general, the correlation functions will 
depend on this scale. However, at a (second-order) critical point the correlation length £ 
diverges, that is the mass vanishes, and there is no length scale in the problem besides 
the overall size L of the domain T>. 

The fact that the scaling limit of (pQ) exists then implies that, instead of starting with a 
lattice model with lattice constant a, we could equally well have started with one with 
some fraction a/b. This would, however, be identical with a lattice model with the original 
spacing a, in which all lengths (including the size of the domain T>) are multiplied by b. 
This implies that the correlation functions in (j2j) are scale covariant 

(0i(6ri)0 2 (6r 2 ) . . .(f)n(br n )) b v = b~^i Xi (^{n)^^) ■ ■ ■ (j> n {r n ))v ■ (5) 

Once again, we can write this in the suggestive form 

^{br) = b-'^ir) , (6) 

as long as what we really mean is (jSj). 

In a massless QFT, the form of the OPE coefficients in (J3J) simplifies: by scale covariance 

where the c^ are pure numbers, and universal if the 2-point functions are normalised so 
that (4>j(ri)(j)j(r2)) = \r\ —r 2 \~ 2xj . (This assumes that the scaling fields are all rotational 
scalars - otherwise it is somewhat more complicated, at least for general dimension.) 

From scale covariance, it is a simple but powerful leap to conformal covariance: suppose 
that the scaling factor b in ((Sj) is a slowly varying function of position r. Then we can try 
to write a generalisation of (jSJ) as 



(Mr[)Mr'2)---Mrn)}v ^J\b{r j )-^{<j> 1 {r 1 )(j> 2 {r 2 )... ( j> n {r n )) 



v 
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where b(r) = \dr'/dr\ is the local jacobian of the transformation r — > r' . 

For what transformations r — > r' do we expect (jSJ) to hold? The heuristic argument runs 
as follows: if the theory is local (that is the interactions in the lattice model are short- 
ranged), then as long as the transformations looks locally like a scale transformation (plus 
a possible rotation), then (jHJ) may be expected to hold. (In Sec. |3]we will make this more 
precise, based on the assumed properties of the stress tensor, and argue that in fact it 
holds only for a special class of scaling fields {0-,} called primary.) 

It is most important that the underlying lattice does not transform (otherwise the state- 
ment is a tautology): (jHJ) relates correlation functions in T>, defined in terms of the limit 
a — > of a model on a regular lattice superimposed on T>, to correlation functions defined 
by a regular lattice superimposed on T>' . 

Transformations which are locally equivalent to a scale transformation and rotation, that 
is, have no local components of shear, also locally preserve angles and are called conformal. 



2.2 Conformal mappings in general 

Consider a general infinitesimal transformation (in flat space) r M — ► r" 1 = r^ + a M (r) 
(we distinguish upper and lower indices in anticipation of using coordinates in which the 
metric is not diagonal.) The shear component is the traceless symmetric part 

all ^d(d + 1) — 1 components of which must vanish for the mapping to be conformal. For 
general d this is very restrictive, and in fact, apart from uniform translations, rotations 
and scale transformations, there is only one other type of solution 

a*{r) = foV - 2(6 • r)^ , 

where b^ is a constant vector. These are in fact the composition of the finite conformal 
mapping of inversion r M — > r M /|r| 2 , followed by an infinitesimal translation 6 M , followed by 
a further inversion. They are called the special conformal transformations, and together 
with the others, they generate a group isomorphic to SO(d + 1, 1). 

These special conformal transformations have enough freedom to fix the form of the 3- 
point functions in R (just as scale invariance and rotational invariance fixes the 2-point 
functions): for scalar operatorgj 

(M r l)M r 2)fo(r 3 )) = |x!+x,-xi I \x*+x,-x, ■ ( 9 ) 

Comparing with the OPE (H1T71) . and assuming non-degeneracy of the scaling dimensions^, 
we see that C123 is the same as the OPE coefficient defined earlier. This shows that the 
OPE coefficients c^ are symmetric in their indices. 

lr The easiest way to show this it to make an inversion with an origin very close to one of the points, 
say n, and then use the OPE, since its image is then very far from those of the other two points. 
2 This and other properties fail in so-called logarithmic CFTs. 
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In two dimensions, the condition that a M (r) be conformal imposes only two differential 
conditions on two functions, and there is a much wider class of solutions. These are more 
easily seen using complex coordinated z = r 1 + ir 2 , z = r 1 — ir 2 , so that the line element 
is ds 2 = dzdz, and the metric is 

a -f° h <r-(° 2 

In this basis, two of the conditions are satisfied identically and the others become 

a z,z = a z,z = ; 

which means that da z /dz = da z /dz = 0, that is, a z is a holomorphic function oc(z) of z, 
and a z is an antiholomorphic function. 

Generalising this to a finite transformation, it means that conformal mappings r — ► r' 
correspond to functions z — > z' = f(z) which are analytic in T>. (Note that the only 
such functions on the whole Riemann sphere are the Mobius transformations f(z) = 
(az + b)/(cz + d), which are the finite special conformal mappings.) 

In passing, let us note that complex coordinates give us a nice way of discussing non- 
scalar fields: if, for example, under a rotation z — > ze ie , (f>j(z, z) -^ e ls i 6 (j)j, we say that (fij 
has conformal spin Sj (not related to quantum mechanical spin), and under a combined 
transformation z — > Xz where A = be" 10 we can write (in the same spirit as (151)) 

^(Xz, Xz) = X-^X^^iz, z) , 

where Xj = Aj + Aj, Sj = Aj — Aj. (Aj, Aj) are called the complex scaling dimensions of 
4>j (although they are usually both real, and not necessarily complex conjugates of each 
other.) 



3 The role of the stress tensor 

Since we wish to explore the consequences of conformal invariance for correlation func- 
tions in a fixed domain T> (usually the entire complex plane), it is necessary to consider 
transformations which are not conformal everywhere. This brings in the stress tensor T^ v 
(also known as the stress-energy tensor or the (improved) energy- momentum tensor). It 
is the object appearing on the right hand side of Einstein's equations in curved space. In 
a classical field theory, it is defined in terns of the response of the action S to a general 
infinitesimal transformation a^(r)\ 

5S=~JT tlv a^d\ (10) 



3 For many CFT computations we may treat z and z as independent, imposing only at the end that 
they should be complex conjugates. 
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(the (l/27r) avoids awkward such factors later on.) Invariance of the action under trans- 
lations and rotations implies that T^ v is conserved and symmetric. Moreover if S is scale 
invariant, T^ v is also traceless. In complex coordinates, the first two conditions imply that 
T z g + T Zz = and T zz = T 2z , so they both vanish, and the conservation equations then 
read d z T zz = 2dT zz /dz = and dT zz /dz = 0. Thus the non-zero components T = T zz 
and T = T zz are respectively holomorphic and antiholomorphic fields. Now if we consider 
a more general transformation for which ot^ u is symmetric and traceless, that is a confor- 
mal transformation, we see that SS = in this case also. Thus, at least classically, we see 
that scale invariance and rotational invariance imply conformal invariance of the action, 
at least if ([TO]) holds. However if the theory contains long-range interactions, for example, 
this is no longer the case. 

In a quantum 2d CFT, it is assumed that the above analyticity properties continue to 
hold at the level of correlation functions: those of T(z) and T(z) are holomorphic and 
antiholomorphic functions of z respectively (except at coincident points.) 

3.0.1 An example - free (gaussian) scalar field 

The prototype CFT is the free, or gaussian, massless scalar field h(r) (we use this notation 
for reasons that will emerge later). It will turn out that many other CFTs are basically 
variants of this. The classical action is 



S[h} = {g/^)j{d tx h){d»h)d 2 r. 

Since h(r) can take any real value, we could rescale it to eliminate the coefficient in front, 
but in later extensions this will have a meaning, so we keep it. In complex coordinates, 
S oc J(d z h)(d z h)d 2 z, and it is easy to see that this is conformally invariant under z — > 
z' = f(z), since d z = f'(z)d z i, d z = f'(z)d z i and d 2 z = \f'(z)\~ 2 d 2 z'. This is confirmed by 
calculating T^ u explicitly: we find T zz = T zz = 0, and 

T = T ZZ = -g(d z h) 2 , T = T- z - z = -g(d- z h) 2 . 

These are holomorphic (resp. antiholomorphic) by virtue of the classical equation of mo- 
tion d z d z h = 0. 

In the quantum field theory, a given configuration {h} is weighted by exp(—S[h]). The 
2-point function isj 

(h(z, z)h(0, 0)) = j J ^ (05 ~ -(1/20) log(^/L 2 ) . 

This means that (T) is formally divergent. It can be made finite, for example, by point- 
splitting and subtracting off the divergent piece: 

T(z) = -ghm^d z h( Z + \5)d z h(z-l5)-^ ] j . (11) 



4 this is cut-off at small k by the assumed finite size L, but we are here also assuming that the points 
are far from the boundary. 
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Figure 1: We consider an infinitesimal transformation which is conformal within C and the 
identity in the complement in T>. 

This doesn't affect the essential properties of T. 



3.1 Conformal Ward identity 

Consider a general correlation function of scaling fields (<f>i(zi, z^fyiizz, z-i) ■ ■ )v in some 
domain T>. We want to make an infinitesimal conformal transformation z — > z' = z + a(z) 
on the points {zj}, without modifying T>. This can be done by considering a contour C 
which encloses all the points {zj} but which lies wholly within T>, such that the trans- 
formation is conformal within C, and the identity z' — z outside C (Fig. [1]). This gives 
rise to an (infinitesimal) discontinuity on C, and, at least classically to a modification of 
the action S according to (fTUj) . Integrating by parts, we find 5S = (l/27r) J c T tJiU a ,J "n l 'd£, 
where n u is the outward-pointing normal and dt is a line element of C. This is more easily 
expressed in complex coordinates, after some algebra, as 

5S = / a(z)T(z)dz + complex conjugate . 

27U Jc 

This extra factor can then be expanded, to first order in a, out of the weight exp(— S[h] — 
SS) ~ (1 — SS) exp(— £[/&]), and the extra piece SS considered as an insertion into the 
correlation function. This is balanced by the explicit change in the correlation function 
under the conformal transformation: 

5((/) 1 (z 1 ,z 1 )(j)2{z2,Z2)...) = — / a{z)(T(z)(f) 1 (z 1 , z 1 )4> 2 {z2, z 2 ) ■■ -)dz + c.c. (12) 

Z7T2 JC 

Let us first consider the case when a(z) = X(z — zi), that is corresponds to a combined 
rotation and scale transformation. In that case 5<fii = (AiA + AiA)0i, and therefore, 
equating coefficients of A and A, 

/dz 
^(Z - Z 1 )(T(z)(f> 1 (z 1 , Z 1 )(j)2(Z2, Z 2 ) ■ ■ ■)— = A 1 ((f) 1 (z 1 ,Z 1 )4>2(z2,Z2) ■ ■ ■) H . 
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Similarly, if we take a = constant, corresponding to a translation, we have S(pj oc d. 

dz 

{'l\z)(p 1 {z 1 , z 1 )(p 2 {z 2 , z 2 ) ■ ■ ■ ) 
ic 



J 3 ^ ^ZjYji 
SO 



^(T{z)M z u zi)4>2(z 2 , z 2 ) ■ ■ ■)— = 2J ^jd Zj ((f>i( z i, h)M z 2, z 2 ) . . .) . 

Using Cauchy's theorem, these two equations tell us about two of the singular terms in 
the OPE of T(z) with a general scaling field 4>j(. z j, z j)'- 

T{z) ■ (pjizj, z j ) = -..+ ,_\ 2 <Pj(zj, Zj) + -—-d^jizj, %) + •••• (13) 

[ZZ-jj Z Z<j 

Note that only integer powers can occur, because the correlation function is a meromorphic 
function of z. 

Example of the gaussian free field. If we take 0' at (r) = e iqh ^ T ' then 

Q 2 /9 



<0!fVi)<»r 2 )> = exp ( - k 2 ((h( ri ) - h{r 2 )f)) r 



K \ri~r 2 \j 

which means that the renormalised field <p q ~ a~ q ' 2s 0' at has scaling dimension x g = q 2 /2g. 
It is then a nice exercise in Wick's theorem to check that the OPE with the stress tensor 
(TT3"j) holds with Aq = x q /2. (Note that in this case the multiplicative renormalisation of 
<f) q is equivalent to ignoring all Wick contractions between fields h(r) at the same point.) 

Now suppose each <ftj is such that the terms 0((z — Zj)~ 2 ~ n ) with n > 1 in (fT3"|) are absent. 
Since a meromorphic function is determined entirely by its singularities, we then know 
the correlation function (T(z) . . .) exactly: 



" \Z-ZjY Z-Zj 7 

This (as well as a similar equation for an insertion of T) is the conformal Ward identity. 
We derived it assuming that the quantum theory could be defined by a path integral and 
the change in the action 5S follows the classical pattern. For a more general CFT, not 
necessarily 'defined' (however loosely) by a path integral, (fill) is usually assumed as a 
property of T. In fact many basic introductions to CFT use this as a starting point. 

Scaling fields 4>j{zj, Zj) such that the most singular term in their OPE with T(z) is 0((z — 
Zj)~ 2 ) are called primary]^ All the other fields like those appearing in the less singular 
terms in (fT3|) are called descendants. Once one knows the correlation functions of all the 
primaries, those of the rest follow from (TTBl rl 



For correlations of such primary fields, we can now reverse the arguments leading to ([13 
for the case of a general infinitesimal conformal transformation ac(z) and conclude that 

S(<f> 1 (zi > Zi)(t) 2 {z 2 , Z 2 ) ...) = J2(^j a '( z j) + <x(Zj)dzj)(<fa.(Z\, Zl)<f>2{z 2 , Z 2 ) ...) , 



5 If we assume that there is a lower bound to the scaling dimensions, such fields must exist. 
6 Since the scaling dimensions of the descendants differ from those of the corresponding primaries by 
positive integers, they are increasingly irrelevant in the sense of the renormalisation group. 
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which may be integrated up to get the result for a finite conformal mapping z — > z' = f(z): 

(</>l(zi, Z 1 )(j) 2 (z 2 , Z 2 ) ...)© = n f'( Z j) A] f'( Z j) ' (M z 'l> z l)M z 2> 4) • • •)& ■ 

j 

This is just the result we wanted to postulate in (jSJ), but now we see that it can hold only 
for correlation functions of primary fields. 

It is important to realise that T itself is not in general primary. Indeed its OPE with 
itself must take the form_| 

T « ■ T <*> = v^t + (T^p™ + ih^™ -■ (15) 

This is because (taking the expectation value of both sides) the 2-point function (T(z)T(zi)) 
is generally non-zero. Its form is fixed by the fact that A T = 2, At = 0, but, since the 
normalisation of T is fixed by its definition (ITU1) . its coefficient c/2 is fixed. This intro- 
duces the conformal anomaly number c, which is part of the basic data of the CFT, along 
with the scaling dimensions (Aj, Aj) and the OPE coefficients Cy&Jj 

This means that, under an infinitesimal transformation a(z), there is an additional term 
in the transformation law of T: 

6T(z) = 2a'(z)T(z) + a(z)d z T(z) + ±ot"(z) . 

For a finite conformal transformation z — > z' = f(z), this integrates up to 

T(z)=f(z) 2 T(z') + f 2 {z' } z} } (16) 

where the last term is the Schwarzian derivative 

w'"(z)w'(z) - \w"{z) 2 



{w(z),z} 



w'(zy 



The form of the Schwarzian can be seen most easily in the example of a gaussian free field. 
In this case, the point-split terms in (fTTl) transform properly and give rise to the first term 
in ([TBI) , but the fact that the subtraction has to be made separately in the origin frame 
and the transformed one leads to an anomalous term 

( f( z + l5)f'(z-±5) 1 \ 

limg 



6^\2g{f{z + \8)-f{z-\8)f 2g8*J 
which, after some algebra, gives the second term in (fl6l) with c = la 



The 0((z — zi)~ 3 ) term is absent by symmetry under exchange of z and z\. 

8 When the theory is placed in curved background, the trace (Tff) ex. cR, where R is the local scalar 
curvature. 

9 This is a classic example of an anomaly in QFT, which comes about because the regularisation 
procedure does not respect a symmetry. 
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Figure 2: The density matrix is given by the path integral over a space-time region in which 
the rows and columns are labelled by the initial and final values of the fields. 



3.2 An application - entanglement entropy 

Let's take a break from the development of the general theory and discuss how the con- 
formal anomaly number c arises in an interesting (and topical) physical context. Sup- 
pose that we have a massless relativistic quantum field theory in 1+1 dimensions, whose 
imaginary time behaviour therefore corresponds to a euclidean CFT. (There are many 
condensed matter systems whose large distance behaviour is described by such a theory.) 
The system is at zero temperature and therefore in its ground state |0), corresponding to 
a density matrix p = |0)(0|. Suppose an observer A has access to just part of the system, 
for example a segment of length C. inside the rest of the system, of total length L 3> £, 
observed by B. The measurements of A and B are entangled. It can be argued that a 
useful measure of the degree of entanglement is the entropy Sa = — Tr^p^logp^, where 
Pa = Tr B p is the reduced density matrix corresponding to A. 

How can we calculate Sa using CFT methods? The first step is to realise that if we can 
compute Tr//^ for positive integer n, and then analytically continue in n, the derivative 
d/dn\ n= i will give the required quantity. For any QFT, the density matrix at finite inverse 
temperature j3 is given by the path integral over some fundamental set of fields h(x, t) (t 
is imaginary time) 



p({h(x,0)},{h(x,i3)}) 



[dh(x,t)]e 



where the rows and columns of p are labelled by the values of the fields at times and f3 
respectively, the the path integral is over all histories h(x, t) consistent with these initial 
and final values (Fig. [2]). Z is the partition function, obtained by sewing together the 
edges at these two times, that is setting h(x,/3) = h(x,0) and integrating J[dh(x, 0)]. 

We are interested in the partial density matrix pa, which is similarly found by sewing 
together the top and bottom edges for x G B, that is, leaving open a slit along the 
interval A (Fig. [3]). The rows and columns of pa are labelled by the values of the fields on 
the edges of the slit. Trp^ is then obtained by sewing together the edges of n copies of 
this slit cylinder in a cyclic fashion (Fig. H|). This gives an n-sheeted surface with branch 
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Figure 3: The reduced density matrix pa is given by the path integral over a cylinder with a 
slit along the interval A. 




Figure 4: Tr p\ corresponds to sewing together n copies so that the edges are connected cycli- 
cally. 

points, or conical singularities, at the ends of the interval A. If Z n is the partition function 
on this surface, then 

Tr p n A = Z n /Z?. 

Let us consider the case of zero temperature, j3 — > oo, when the whole system is in the 
ground state |0). If the ends of the interval are at (xi,x 2 ), the conformal mapping 

w — X\ \ 1 / n 
,w — x 2 J 

maps the n-sheeted w-surface to the single-sheeted complex z-plane. We can use the 
transformation law ( fl~6|) to compute (T(w)), given that (T(z)) = by translational and 
rotational invariance. This gives, after a little algebra, 

(c/12){l-l/n*){x 2 - Xl )* 

( T M) = 7 77J7 77J • 

(w — X\Y\W — x 2 ) 

Now suppose we change the length £ = \x 2 — Xi\ slightly, by making an infinitesimal 
non-conformal transformation x — > x + 5£5(x — Xq), where x\ < Xq < x 2 . The response of 
the log of the partition function, by the definition of the stress tensor, is 

6\ogZ n = -—— (T xx (x ,t))dt 

(the factor n occurs because it has to be inserted on each of the n sheets.) Writing 
T xx = T + T, the integration in each term can be carried out by wrapping the contour 
around either of the points X\ or x 2 . The result is 

<91ogZ„ (c/6)(w-l/ra) 

dl I 

so that Z n /Z™ ~ £-( c / 6 )( n - 1 / n ). Taking the derivative with respect to n at n = 1 we get 
the final result 

S A ~(c/3)log£. 
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Figure 5: The state \<pj) is defined by weighting field configurations on the circle with the path 
integral inside, with an insertion of the field (ftj(0) at the origin. 

4 Radial quantisation and the Virasoro algebra 

4.1 Radial quantisation 

Like any quantum theory CFT can be formulated in terms of local operators acting on a 
Hilbert space of states. However, as it is massless, the usual quantisation of the theory on 
the infinite line is not so useful since it is hard to disentangle the continuum of eigenstates 
of the hamiltonian, and we cannot define asymptotic states in the usual way. Instead it 
is useful to exploit the scale invariance, rather than the time-translation invariance, and 
quantise on a circle of fixed radius r . In the path integral formulation, heuristically the 
Hilbert space is the space of field configurations \{h(ro, 8)}) on this circle. The analogue 
of the hamiltonian is then the generator D of scale transformations. It will turn out that 
the spectrum of D is discrete. In the vacuum state each configuration is weighted by the 
path integral over the disc \z\ < r , conditioned on taking the assigned values on r , see 
Fig. U 

|0)= [{dh(r<r )]e- s W\{h(r ,6)}). 



The scale invariance of the action means that this state is independent of r , up to a 
constant. If instead we insert a scaling field (fij(0) into the above path integral, we get 
a different state \<f)j)- On the other hand, more general correlation functions of scaling 
fields are given in this operator formalism by 10 l 

((j) j (z 1 ,Z 1 )(j) j (Z2,Z2)) = (O|R0j(Zi, z 1 )(j)j(z 2 , z 2 )\0) , 

where R means the r-ordered product (like the time-ordered product in the usual case), 
and (0| is defined similarly in terms of the path integral over r > r . Thus we can identify 



\<j)j) = hmfaiz, z)\0) . 

z — >U 

This is an example of the operator- state correspondence in CFT. 



10 



We shall make an effort consistently to denote actual operators (as opposed to fields) with a hat. 
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Just as the hamiltonian in ordinary QFT can be written as an integral over the appropriate 
component of the stress tensor, we can write 



JL. 1 f - 

D = L + L = / zT(z)dz + c.c 

2m Jc 



2m' Jc 

where C is any closed contour surrounding the origin. This suggests that we define more 
generally 

L n = ^- I z n+1 f(z)dz, 

2m JC 

and similarly L n . 

If there are no operator insertions inside C it can be shrunk to zero for n > — 1, thus 

L n |0) = forn>-l. 

On the other hand, if there is an operator fy inserted at the origin, we see from the OPE 
flEU that 

Lol^-) = A il&> • 

If 4>j is primary, we further see from the OPE that 

L n \(j)j) =0 for n > 1 , 
while for n < — 1 we get states corresponding to descendants of <pj. 

4.2 The Virasoro algebra 

Now consider L m L n acting on some arbitrary state. In terms of correlation functions this 
involves the contour integrals 

dz2 m+ i f dzi +lr 



where C2 lies outside C\ , because of the R-ordering. If instead we consider the operators 
in the reverse order, the contours will be reversed. However we can then always distort 
them to restore them to their original positions, as long as we leave a piece of the z 2 
contour wrapped around z\. This can be evaluated using the OPE (IT51) of T with itself: 

dZl n+1 l dz 2 +1 / 2 rrf_\ , 1 o rrf„\ , c / 2 \ 



JC 27U J 2m \ [Zo — Z\) Zo — Z] 



\z 2 -z x y z 2 -zi )z 2 -z 1 y) 

= f ^? +1 (2(m + l)*;T(*i) + z^^i^i) + 75™K - 1)^™- 2 ) • 

Jc 2lTl v " / 

The integrals can then be re-expressed in terms of the L n . This gives the Virasoro algebra: 
[L m , L n ] = (m- n)L m+n + ^m(m 2 - l)5 m+nfi , (17) 
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with an identical algebra satisfied by the L n . It should be stressed that ( TT71) is completely 
equivalent to the OPE (TT5l) . but of course algebraic methods are often more efficient in 
understanding the structure of QFT. The first term on the right hand side could have been 
foreseen if we think of L n as being the generator of infinitesimal conformal transformations 
with a(z) oc z n+1 . Acting on functions of z this can therefore be represented by £ n = 
z n+1 d z , and it is easy to check that these satisfy (JT7J) without the central term (called the 
Witt algebra.) However the L n act on states of the CFT rather than functions, which 
allows for the existence of the second term, the central term. The form of this, apart from 
the undetermined constant c, is in fact dictated by consistency with the Jacobi identity. 
Note that there is a closed subalgebra generated by (I/ 1; L , L-i), which corresponds to 
special conformal transformations. 

One consequence of (IT7|) is 

[Lq, L_ n ] = nL_ n , 

so that Lo(L- n \(j)j)) = (Aj+n)L^ n \<f)j) , which means that the L n with n < act as raising 
operators for the weight, or scaling dimension, A, and those with n > act as lowering 
operators. The state \<f>j) corresponding to a primary operator is annihilated by all the 
lowering operators. It is therefore a lowest weight statel^i By acting with all possible 
raising operators we build up a lowest weight representation (called a Verma module) of 
the Virasoro algebra: 



£-3|0i),£-2£-i|0j),i£il0i); 

L-al^Llilfc); 

£-i|&>; 

\4>j) • 

4.3 Null states and the Kac formula 

One of the most important issues in CFT is whether, for a given c and Aj, this repre- 
sentation is unitary, and whether it is reducible (more generally, decomposable). It turns 
out that these two are linked, as we shall see later. Decomposability implies the existence 
of null states in the Verma module, that is, some linear combination of states at a given 
level is itself a lowest state. The simplest example occurs at level 2, if 



L n (L_ 2 \<f )j )-(l/g)L 2 _ 1 \<P J ))=0. 



for n > (the notation with g is chosen to correspond to the Coulomb gas later on.) By 
taking n — 1 and n = 2 and using the Virasoro algebra and the fact that \<f>j) is a lowest 



ii 



In the literature this is often called, confusingly, a highest weight state. 
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weight state, we get 

A _ 3£~ 2 _ (3g-2)(3-2g) 

4 5 

This is special case (r, s) = (2, 1) of the Kac formula: with c parametrised as above, i i 12 l 



A, = A r , s (g) = 



±g 



then \<f>j) has a null state at level r • s. We will not prove this here, but later will indicate 
how it is derived from Coulomb gas methods. 

Removing all the null states from a Verma module gives an irreducible representation of 
the Virasoro algebra. Null states (and all their descendants) can consistently be set to 
zero in a given CFT. (This is no guarantee that they are in fact absent, however.) 

One important consequence of the null state is that the correlation functions of <pj{z, z) 
satisfy linear differential equations in z (or z) of order rs. The case rs = 2 will be 
discussed as an example in the last lecture. This allows us in principle to calculate all the 
four-point functions and hence the OPE coefficients. 



4.4 Fusion rules 

Let us consider the 3-point function 

where the first two fields sit at the indicated places in the Kac table, but the third is a 
general primary scaling field of dimension A. The form of this is given by ([2]). If we insert 
J c (z — Zi)~ 1 T(z)dz into this correlation function, where C surrounds z\ but not the other 
two points, this projects out L_ 2 4 ) 2,i oc ^02,1- On the other hand, the full expression is 
given by the Ward identity ffT4|) . After some algebra, we find that this is consistent only 
if 

A = A r± i iS , 

otherwise the 3-point function, and hence the OPE coefficient of a in 02, i • <Pr,s, vanishes. 

This is an example of the fusion rules in action. It shows that Kac operators compose 
very much like irreducible representations of SU(2), with the r label playing the role 
of spin |(r — 1). The s-labels compose in the same way. More generally the fusion rule 
coefficients Nf, tell us not only which OPEs can vanish, but which ones actually do appear 
in a particular CFT£j In this simplest case we have (suppressing the s-indices for clarity) 

N^ r i = <5 r » iT . +r '_i + 5 r .» iT . +r '_3 + • • ■ + Or",!? — T-'j-fl • 



12 Note that g and 1/g given the same value of c, and that A rtS (g) = A Str (l/g). This has led to endless 
confusion in the literature. 

13 They can actually take values > 2 if there are distinct primary fields with the same dimension. 
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A very important thing happens if g is rational = p/p'. Then we can write the Kac 

formula as 

(rp — sp') 2 — (p — p') 2 

4pp' 
and we see that A rjS = Ap/_ riP _ s , that is, the same primary field sits at two different 
places in the rectangle 1 < r < p' — 1, 1 < s < p — 1, called the Kac table. If we now 
apply the fusion rules to these fields we see that we get consistency between the different 
constraints only if the fusion algebra is truncated, that fields within the rectangle do not 
couple to those outside. 

This shows that, at least at the level of fusion, we can have CFTs with a finite number 
of primary fields. These are called the minimal models 1 !. However, it can be shown that, 
among these, the only ones admitting unitary representations of the Virasoro algebra, that 
is for which (ijj\if)) > for all states \ip) in the representation, are those with \p — p'\ = 1 
and p,p' > 3. Moreover these are the only unitary CFTs with c < 1. The physical 
significance of unitarity will be mentioned shortly. 



5 CFT on the cylinder and torus 

5.1 CFT on the cylinder 

One of the most important conformal mappings is the logarithmic transformation w = 
(L/2ir) log z, which maps the z-plane (punctured at the origin) to an infinitely long cylin- 
der of circumference L (equivalently, a strip with periodic boundary conditions.) It is 
useful to write w = t + iu, and to think of the coordinate t running along the cylinder as 
imaginary time, and u as space. CFT on the cylinder then corresponds to euclidean QFT 
on a circle. 

The relation between the stress tensor on the cylinder and in the plane is given by ( flTfl) : 

T(w) cyl = (dz/dw) 2 T(z) + ^{z,w} = (2tt/L) 2 (z 2 T(z) plane - ^) , 

where the last term comes from the Schwarzian derivative. 

The hamiltonian H on the cylinder, which generates infinitesimal translations in t, can be 
written in the usual way as an integral over the time-time component of the stress tensor 

H = ±- f L f tt (u)du = ±- f L (f(u)+f(u))du, 

Z7T JO Z71 JO 

which corresponds in the plane to 



14 The minimal models are examples of rational CFTs: those which have only a finite number of fields 
which are primary with respect to some algebra, more generally an extended one containing Virasoro as 
a sub-algebra. 



CFT and Statistical Mechanics 19 

Similarly the total momentum P, which generates infinitesimal translations in u, is the 
integral of the T tu component of the stress tensor, which can be written as (2n/L)(L — Lq). 

Eq. (I19p . although elementary, is one of the most important results of CFT in two dimen- 
sions. It relates the dimensions of all the scaling fields in the theory (which, recall, are the 

eigenvalues of Lq and Lq) to the spectra of H and P on the cylinder. If we have a lattice 
model on the cylinder whose scaling limit is described by a given CFT, we can therefore 
read off the scaling dimensions, up to finite-size corrections in (a/L), by diagonalising the 
transfer matrix t ~ 1 — aH. This can be done either numerically for small values of L, 
or, for integrable models, by solving the Bethe ansatz equations. 

In particular, we see that the lowest eigenvalue of H (corresponding to the largest eigen- 
value of the transfer matrix) is 

„ tic 2n , „ - r . 

E = -— + —(A + A ), 

where (Ao, A ) are the lowest possible scaling dimensions. In many CFTs, and all unitary 
ones, this corresponds to the identity field, so that A = A = 0. This shows that c can 
be measured from finite-size behaviour of the ground state energy. 

Eq also gives the leading term in the partition function Z = Tr e~ eH on a finite cylinder 
(a torus) of length £ 3> L. Equivalently, the free energy (in units of ksT) is 

,-, -, ^ kcE 

F = -logZ~-—. 
bL 

In this equation F represents the scaling limit of free energy of a 2d classical lattice 
modeO However we can equally well think of t as being space and u imaginary time, in 
which case periodic boundary conditions imply finite inverse temperature j3 = 1/kgT = L 
in a Id quantum field theory. For such a theory we then predict that 

ixclksL 



6 ' 

or, equivalently, that the low-temperature specific heat, at a quantum critical point de- 
scribed by a CFT (generally, with a linear dispersion relation uo ~ \q\), has the form 

nck 2 B T 
G,~-3-. 

Note that the Virasoro generators can be written in terms of the stress tensor on the 
cylinder as 

L n = ^- [ L e inu f{u,0)du. 

Ztt Jo 



15 This treatment overlooks UV divergent terms in F of order (£L/a 2 ), which are implicitly set to zero 
by the regularisation of the stress tensor. 
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1 + 1 




1 

Figure 6: A general torus is obtained by identifying opposite sides of a parallelogram. 

In a unitary theory, T is self-adjoint, and hence U n = L_ n . Unitarity of the QFT corre- 
sponds to reflection positivity of correlation functions: in general 

(0i (ui, ti)<j) 2 (u 2 , t 2 )... (j>i (iti, -£i)0 2 (M2, -t 2 )) 

is positive if the transfer matrix T can be made self-adjoint, which is generally true if the 
Boltzmann weights are positive. Note however, that a given lattice model (e.g. the Ising 
model) contains fields which are the scaling limit of lattice quantities in which the transfer 
matrix can be locally expressed (e.g. the local magnetisation and energy density) and for 
which one would expect reflection positivity to hold, and other scaling fields (e.g. the 
probability that a given edge lies on a cluster boundary) which are not locally expressible. 
Within such a CFT, then, one would expect to find a unitary sector - in fact in the Ising 
model this corresponds to the p = 3,p' = 4 minimal model - but also possible non-unitary 
sectors in addition. 

5.2 Modular invariance on the torus 

We have seen that unitarity (for c < 1), and, more generally, rationality, fix which scaling 
fields may appear in a given CFT, but they don't fix which ones actually appear. This 
is answered by considering another physical requirement: that of modular invariance on 
the torus. 

We can make a general torus by imposing periodic boundary conditions on a parallelo- 
gram, whose vertices lie in the complex plane. Scale invariance allows us to fix the length 
of one of the sides to be unity: thus we can choose the vertices to be at (0, 1, 1 + r, r), 
where r is a complex number with positive imaginary part. In terms of the conventions 
of the previous section, we start with a finite cylinder of circumference L = 1 and length 
Imr, twist one end by an amount Rer, and sew the ends together - see Fig. An 
important feature of this parametrisation of the torus is that it is not unique: the trans- 
formations T : t —* t + 1 and S : r — > — 1/r give the same torus (Fig. [7j). Note that 
S interchanges space u and imaginary time t in the QFT. S and T generate an infinite 

discrete group of transformations 

ar + b 

ct + a 
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Figure 7: Two ways of viewing the same torus, corresponding to the modular transformations 
T and S. 

with (a, b, c, d) all integers and ad — be = 1. This is called SL(2, Z) or the modular group. 
Note that S 2 = 1 and (ST) 3 = 1. 

Consider the scaling limit of the partition function Z of a lattice model on this torus. 
Apart from the divergent term in logZ, proportional to the area divided by a 2 , which in 
CFT is set to zero by regularisation, the rest should depend only on the aspect ratio of 
the torus and thus be modular invariant. This would be an empty statement were it not 
that Z can be expressed in terms of the spectrum of scaling dimensions of the CFT in a 
manner which is not, manifestly, modular invariant. 

Recall that the generators of infinitesimal translations along and around the cylinder can 
be written as 

H = 2tt(L + 1 ) - f P = 2n(L - I ) . 

The action of twisting the cylinder corresponds to a finite translation around its circum- 
ference, and sewing the ends together corresponds to taking the trace. Thus 

Z = r x re -( Im ^+*(R e - r )'P 

_ TTclmr/6 rp 2-KirLo -2-rriLo 

= (qq)- c/24 Trq Lo q To , 
where in the last line we have defined q = e 2mT . 

The trace means that we sum over all eigenvalues of L$ and Lq, that is all scaling fields 
of the CFT. We know that these can be organised into irreducible representations of 
the Virasoro algebra, and therefore have the form (A + N, A + N), where A and A 
correspond to primary fields and (TV, TV) are non-negative integers labelling the levels of 
the descendants. Thus we can write 

Z = Y L n A - K XA{q)XA^) . 

A, A 

where n AA is the number of primary fields with lowest weights (A, A), and 



XA(q) = g" c/24+A E ^(TV)g 

N=0 



N 
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where cIa(N) is the degeneracy of the representation at level N. It is purely a property 
of the representation, not the particular CFT, and therefore so is Xa(q)- This is called 
the character of the representation. 

The requirement of modular invariance of Z under T is rather trivial: it says that all fields 
must have integer conformal spin 16 !. However the invariance under S is highly non-trivial: 
it states that Z, which is a power series in q and q, can equally well be expressed as an 
identical power series in q = e~ 2m l T and q. 

We can get some idea of the power of this requirement by considering the limit q — *■ 1, 
q — > 0, with q real. Suppose the density of scaling fields (including descendants) with 
dimension x = A + A in the range (x, x + 5x) (where 1 ^> 5x ^> x) is p(x)5x. Then, in 
this limit, when q — 1 — e, e -C 1, 



/oo /*oo 

p(x)e xlogq dx ~ p(x)e- £X dx. 



On the other hand, we know that as q — > 0, Z ~ g c/i2+x ^ e (27r) (c/12 z )A w here x < 
is the lowest scaling dimension (usually 0). Taking the inverse Laplace transform, 

p(x) ~ /" e ^+(2T) 2 (c/12-x )/6_^_ 

Using the method of steepest descents we then see that, as x — > oo, 

p{x) ~ exp (4tt(c/12 - x ) 1/2 x 1/2 ) , 

times a (calculable) prefactor. This relation is of importance in understanding black hole 
entropy in string theory. 



5.2.1 Modular invariance for the minimal models 

Let us apply this to the minimal models, where there is a finite number of primary 
fields, labelled by (r, s). We need the characters Xr,s(o)- If there were no null states, the 
degeneracy at level iV would be the number of states of the form . . . L^L^L^W^) with 
J2j J n j — N. This is just the number of distinct partitions of iV into positive integers, 
and the generating function is rifcli(l — q k )~ l ■ 

However, we know that the representation has a null state at level rs, and this, and all 
its descendants, should be subtracted off. Thus 



J ; 



xrs(q)= q - c/M m i - ( i k )- i a-q rs + 



fc=l 



But, as can be seen from the Kac formula (j!8j) . A rjS + rs = A p / +r . )P _ s , and therefore the 
null state at level rs has itself null states in its Verma module, which should not have been 



16 It is interesting to impose other kinds of boundary conditions, e.g. antiperiodic, on the torus, when 
other values of the spin can occur. 
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subtracted off. Thus we must add these back in. However, it is slightly more complicated 
than this, because for a minimal model each primary field sits in two places in the Kac 
rectangle, A rjS = A p '_ r>p _ s . Therefore this primary field also has a null state at level 
[p 1 — r)(p — s), and this has the dimension A 2p '- r ,s = A ri2p -s and should therefore also 
be added back in if it has not be included already. A full analysis requires understanding 
how the various submodules sit inside each other, but fortunately the final result has a 
nice form 

oo 

Xrs(q) = q- c/24 m i ~ V k T l ( K rM) - Kr t -M)) , W 

fe=l 

where 

oo 

K r>a {q) = Y, q( 2npp ' +rp - sp '^/ 4pp ' . (21) 

n=— oo 

The partition function can then be written as finite sum 

" / j n r,s;f,sXrs[Q)Xfs[C[) / , n r,s;f,sXrs\Q)Xrs\C[) ■ 

r,s;f,s r,s;r,s 

The reason this can happen is that the characters themselves transform linearly under 
S : q — > q, as can be seen (after quite a bit of algebra, by applying the Poisson sum 
formula to (J21D and Euler's identities to the infinite product): 

Xrs(q) = J2 S rs S 'Xr's'(q), 
r',s' 

where S is a matrix whose rows and columns are labelled by (rs) and (rV). Another way 
to state this is that the characters form a representation of the modular group. The form 
of S is not that important, but we give it anyway: 



1/2 , _ , _ 

'rs 



s:'z = — (_iy+rs + sf sin >2^j_ sin 



y pp' J p' p 

The important properties of S are that it is real and symmetric and S 2 = 1. 

This immediately implies that the diagonal combination, with n rjS . f)5 = <5 r f<5 s5 , is modular 
invariant: 

Y,Xrs(q)Xrs(q) =J2J2J2 S r r ' s S 'S 7 r 's S "x r > s >(q)Xr"s"(q) = J2 Xrs(q)Xrs(q) , 
r,s r,s r ',s' r",s" r,s 

where we have used SS T = 1. This gives the diagonal series of CFTs, in which all possible 
scalar primary fields in the Kac rectangle occur just once. These are known as the A n 
series. 

It is possible to find other modular invariants by exploiting symmetries of S. For example, 
if p'/2 is odd, the space spanned by Xrs + Xp'-r,s, with r odd, is an invariant subspace 
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of S, and is multiplied only by a pure phase under T. Hence the diagonal combination 
within this subspace 

•^ / y \Xr,s + Xp'—r,s\ 

r odd,s 

is modular invariant. Similar invariants can be constructed if p/2 is odd. This gives the 
D n series. Note that in this case some fields appear with degeneracy 2. Apart from these 
two infinite series, there are three special values of p and p' (12,18,30) denoted by E 6j 7 i8 . 
The reason for this classification will become more obvious in the next section when we 
construct explicit lattice models which have these CFTs as their scaling limits] 17 ! 



6 Height models, loop models and Coulomb gas meth- 
ods 

6.1 Height models and loop models 

Although the ADE classification of minimal CFTs with c < 1 through modular invariance 
was a great step forwards, one can ask whether there are in fact lattice models which have 
these CFTs as their scaling limit. The answer is yes - in the form of the ADE lattice 
models. These can be analysed non-rigorously by so-called Coulomb gas methods. 

For simplicity we shall describe only the so-called dilute models, defined on a triangular 
latticeo At each site r of the lattice is defined a 'height' h{r) which takes values on 
the nodes of some connected graph Q. An example is the linear graph called A m shown 
in Fig. [HI in which h(r) can be thought of as an integer between 1 and m. There is a 



Figure 8: The graph A m , with m nodes. 

restriction in these models that the heights on neighbouring sites of the triangular lattice 
must either be the same, or be adjacent on Q. It is then easy to see that around a given 
triangle either all three heights are the same (which carries relative weight 1), or two 
of them are the same and the other is adjacent on Q^j In this case, if the heights are 
(h, h', h'), the weight is x(Sh/Sh'Y^ 6 where Sh is a function of the height h, to be made 
explicit later, and x is a positive temperature- like parameter o (A simple example is A2, 



17 This classification also arises in the finite subgroups of SU(2), of simply-laced Lie algebras, and in 
catastrophe theory. 

18 Similar models can be defined on the square lattice. They give rise to critical loop models in the 
dense phase. 

19 Apart from the pathological case when Q itself has a 3-cycle, in which case we can enforce the 
restriction by hand. 

20 It is sometimes useful, e.g. for implementing a transfer matrix, to redistribute these weights around 
the loops. 
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Figure 9: If the heights on the vertices are not all equal, we denote this by a segment of a curve 
on the dual lattice, as shown. 

corresponding to the Ising model on the triangular lattice.) 

The weight for a given configuration of the whole lattice is the product of the weights 
for each elementary triangle. Note that this model is local and has positive weights if Sh 
is a positive function of h. Its scaling limit at the critical point should correspond to a 
unitary CFT. 

The height model can be mapped to a loop model as follows: every time the heights in a 
given triangle are not all equal, we draw a segment of a curve through it, as shown in Fig. [91 
These segments all link up, and if we demand that all the heights on the boundary are 
the same, they form a set of nested, non-intersecting closed loops on the dual honeycomb 
lattice, separating regions of constant height on the original lattice. Consider a loop for 
which the heights just inside and outside are h and h! respectively. The loop has convex 
(outward-pointing) and concave (inward-pointing) corners. Each convex corner carries a 
factor (Sh/Sh') 1 ^ 6 , and each concave corner the inverse factor. But each loop has exactly 
6 more outward pointing corners than inward pointing ones, so it always carries an overall 
weight Sh/S' h , times a factor x raised to the length of the loop. Let us now sum over the 
heights consistent with a fixed loop configuration, starting with innermost regions. Each 
sum has the form 

2J {s h /s h i) , 

h:\h-h'\=l 

where \h—h'\ = 1 means that h and h! are adjacent on Q. The next stage in the summation 
will be simple only if this is a constant independent of h! . Thus we assume that the Sh 
satisfy 

2 , Sh = ASh> , 

h:\h-h'\=l 

that is, Sh is an eigenvector of the adjacency matrix of Q, with eigenvalue A. For A m , for 

example, these have the form 

nkh 
S h oc sm — , (22) 

771 + 1 

where 1 < k < m, corresponding to A = 2 cos (irk/(m + 1)). Note that only the case 
k — 1 gives all real positive weights. 
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Figure 10: Nested set of loops, separating regions of constant height. We iteratively sum over 
the heights h, starting with the innermost. 



Having chosen the Sh in this way, we can sum out all the heights consistent with a 
given loop configuration (Fig. [TU]) . starting with the innermost and moving outwards, and 
thereby express the partition function as a sum over loop configurations: 

17 _ \ " a number of loops total length (O'i') 

loop configs 

When x is small, the heights are nearly all equal (depending on the boundary condition,) 
and the typical loop length and number is small. At a critical point x = x c we expect 
this to diverge. Beyond this, we enter the dense phase, which is still critical in the loop 
sense, even though observables which are local in the original height variables may have 
a finite correlation length. For example, for x > x c in the Ising model, the Ising spins are 
disordered but the the cluster boundaries are the same, in the scaling limit, as those of 
critical percolation for site percolation on the triangular lattice. 

However, we could have obtained the same expression for Z in several different ways. 
One is by introducing n-component spins s a (R) with a = 1, . . . , n on the sites of the dual 
lattice, and the partition function 



'0(n) 



Trn 

RR' 



X 



Y,S a (R)s a (R') 



o=l 



where the product is over edges of the honeycomb lattice, the trace of an odd power of 
s a (R) is zero, and Trs a (R)sb(R) = S a b. Expanding in powers of x, and drawing in a curve 
segment each time the term proportional to x is chosen on a given edge, we get the same 
set of nested loop configurations, weighted as above with A = n. This is the O(n) model. 
Note that the final expression makes sense for all real positive values of n, but it can 
be expressed in terms of weights local in the original spins only for positive integer n. 
Only in the latter case do we therefore expect that correlations of the O(n) spins will 
satisfy reflection positivity and therefore correspond to a unitary CFT, even though the 
description in terms of heights is unitary. This shows how different sectors of the 'same' 
CFT can describe rather different physics. 
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Figure 11: Non-contractible loops on the cylinder can be taken into account by the insertion of 
a suitable factor along a seam. 

6.2 Coulomb gas methods 

These loop models can be solved in various ways, for example by realising that their 
transfer matrix gives a representation of the Temperley-Lieb algebra, but a more powerful 
if less rigorous method is to use the so-called Coulomb gas approach. Recalling the 
arguments of the previous section, we see that yet another way of getting to (123|) is by 
starting from a different height model, where now the heights h(r) are defined on the 
integers (times ir, for historical reasons) 7rZ. As long as we choose the correct eigenvalue 
of the adjacency matrix, this will give the same loop gas weights. That is we take Q to 
be Aoo. In this case the eigenvectors correspond to plane wave modes propagating along 
the graph, labelled by a quasi- momentum \ with |%| < 1: Sh oc e txh , whose eigenvalue is 
A = 2cos(7rx)- Because these modes are chiral, we have to orient the loops to distinguish 
between \ an d — X- Each oriented loop then gets weighted with a factor e ±J7rx/ ' 6 at each 
vertex of the honeycomb lattice it goes through, depending on whether it turns to the left 
or the right. 

This version of the model, where the heights are unbounded, is much easier to analyse, 
at least non-rigorously. In particular, we might expect that in the scaling limit, after 
coarse-graining, we can treat hir) as taking all real values, and write down an effective 
field theory. This should have the property that it is local, invariant under h(r) — > h(r) + 
constant, and with no terms irrelevant under the RG (that is, entering the effective action 
with positive powers of a.) The only possibility is a free gaussian field theory, with action 



S 



(g /4n) J(Vh) 2 d 2 r . 



However, this cannot be the full answer, because we know this corresponds to a CFT with 
c = 1. The resolution of this is most easily understood by considering the theory on a long 
cylinder of length £ and circumference L <ti £. Non-contractible loops which go around 
the cylinder have the same number of inside and outside corners, so they are incorrectly 
counted. This can be corrected by inserting a factor fit e txh ^' ' e~ lxh ^ t+1 '°\ which counts 
each loop passing between (£, 0) and (t + 1, 0) with just the right factors e ±mx . These 
factors accumulate to e* x/l< - _ ^ //2,0 - ) e~ %xh ^l 2,Q \ corresponding to charges ±x at the ends of 
the cylinder. This means that the partition function is 

Z ~ Z c=1 (e ixh ^e- ixh( - e ^) . 
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But we know (Sec. 13.11) that this correlation function decays like r~ 2Xx in the plane, where 
x q = q 2 /2g , and therefore on the cylinder 

Z ^ e Tre/6L e -27T( x 2 /2g )£/L 

from which we see that the central charge is actually 

6X 2 



c= 1 



90 



However, we haven't yet determined go- This is fixed by the requirement that the screening 
fields e ±l2h( - r ' , which come from the the fact that originally h(r) G 7rZ, should be marginal, 
that is they do not affect the scaling behaviour so that we can add them to the action 
with impunity. This requires that they have scaling dimension X2 = 2. However, now x q 
should be calculated from the cylinder with the charges e ±l xM±^/ 2 ) a ^ the ends: 

= (g±x) 2 _ x^ = q 2 ± 2 xq / 24 x 

9 go go go 

Setting X2 = 2 we then find go = 1 ± % and therefore 

c-l- 6 '^- 1 ' 2 . 

go 

6.3 Identification with minimal models 

The partition function for the height models (at least on cylinder) depends only on the 

eigenvalue A of the adjacency matrix and hence the Coulomb gas should work equally 

well for the models on A m if we set x = k/(m + 1). The corresponding central charge is 

then 

Qk 2 
c = 1 — 



(m + l)(m+ l±k) ' 
If we compare this with the formula for the minimal models 

6(p — p') 2 



1- 



pp' 



we are tempted to identify k = p — p' and m + 1 = p'. This implies go = v/ft \ which can 
therefore be identified with the parameter g introduced in the Kac formula] 21 ! Moreover, 
if we compute the scaling dimensions of local fields <f> r (R) = cos ((r — l)kh(R)/(m + IV) 
using ff24l) we find perfect agreement with the leading diagonal A rr . of the Kac table] 22 ! 
We therefore have strong circumstantial evidence that the scaling limit of the dilute A p >_i 



21 The other solution corresponds to interchanging p' and p. 
22 These are the relevant fields in the RG sense. 
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models (choosing the eigenvalue A = 2cos(ir(p — p')/p')) is the (p,p') minimal model with 
p > p'. Note that only if k = 1, that is p = p' + 1, are these CFTs unitary, and this is 
precisely the case where the weights of the lattice model are real and positive. 

For other graphs Q we can try to make a similar identification. However, this is going 
to work only if the maximal eigenvalue of the adjacency matrix of Q is strictly less than 
2. A famous classification then shows that this restricts Q to be either of the form A m , 
D m or one of three exceptional cases £"6,7,80 These other graphs also have eigenvalues of 
the same form (refeigen) as A m , but with m + 1 now being the Coxeter number, and the 
allowed integers k being only a subset of those appearing in the A-series. These correspond 
to the Kac labels of the allowed scalar operators which appear in the appropriate modular 
invariant partition function. 



7 Boundary conformal field theory 

7.1 Conformal boundary conditions and Ward identities 

So far we haven't considered what happens at the boundary of the domain T>. This is a 
subject with several important applications, for example to quantum impurity problems 
(see the lectures by Affleck) and to D-branes in string theory. 

In any field theory in a domain with a boundary, one needs to consider how to impose 
a set of consistent boundary conditions. Since CFT is formulated independently of a 
particular set of fundamental fields and a lagrangian, this must be done in a more general 
manner. A natural requirement is that the off-diagonal component Tj|j_ of the stress 
tensor parallel/perpendicular to the boundary should vanish. This is called the conformal 
boundary condition. If the boundary is parallel to the time axis, it implies that there is 
no momentum flow across the boundary. Moreover, it can be argued that, under the RG, 
any uniform boundary condition will flow into a conformally invariant one. For a given 
bulk CFT, however, there may be many possible distinct such boundary conditions, and 
it is one task of BCFT to classify these. 

To begin with, take the domain to be the upper half plane, so that the boundary is the 
real axis. The conformal boundary condition then implies that T(z) = T(z) when z is on 
the real axis. This has the immediate consequence that correlators of T are those of T 
analytically continued into the lower half plane. The conformal Ward identity now reads 



23 A = 2 corresponds to the extended diagrams A m , etc., which give interesting rational CFTs with 
c = 1. However models based on graphs with A > 2 probably have a different kind of a transition at 
which the mean loop length remains finite. 
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■j K Zj,Zj)). (25) 

In radial quantisation, in order that the Hilbert spaces defined on different hypersurfaces 
be equivalent, one must now choose semicircles centered on some point on the boundary, 
conventionally the origin. The dilatation operator is now 

b =hL zf w z -hL>*w> (26) 

where S is a semicircle. Using the conformal boundary condition, this can also be written 
as 

D = L = -^ / zf{z)dz , (27) 

2m Jc 

where C is a complete circle around the origin. 

Note that there is now only one Virasoro algebra. This is related to the fact that con- 
formal mappings which preserve the real axis correspond to real analytic functions. The 
eigenstates of L correspond to boundary operators <f>j(0) acting on the vacuum state |0). 
It is well known that in a renormalizable QFT fields at the boundary require a different 
renormalization from those in the bulk, and this will in general lead to a different set of 
conformal weights. It is one of the tasks of BCFT to determine these, for a given allowed 
boundary condition. 

However, there is one feature unique to boundary CFT in two dimensions. Radial quan- 
tization also makes sense, leading to the same form (12 7p for the dilation operator, if the 
boundary conditions on the negative and positive real axes are different. As far as the 
structure of BCFT goes, correlation functions with this mixed boundary condition be- 
have as though a local scaling field were inserted at the origin. This has led to the term 
'boundary condition changing (bcc) operator'. 

7.2 CFT on the annulus and classification of boundary states 

Just as consideration of the partition function on the torus illuminates the bulk operator 
content n A A , it turns out that consistency on the annulus helps classify both the allowed 
boundary conditions, and the boundary operator content. To this end, consider a CFT 
in an annulus formed of a rectangle of unit width and height 5, with the top and bottom 
edges identified (see Fig. [12]). The boundary conditions on the left and right edges, labelled 
by a, b, . . ., may be different. The partition function with boundary conditions a and b on 
either edge is denoted by Z a b(S). 

One way to compute this is by first considering the CFT on an infinitely long strip of unit 
width. This is conformally related to the upper half plane (with an insertion of boundary 
condition changing operators at and oo if a ^ b) by the mapping z — > (1/V) In z. The 
generator of infinitesimal translations along the strip is 

H ab = txD- nc/24 = tvLq - vrc/24 . (28) 
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Figure 12: The annulus, with boundary conditions a and b on either boundary. 



Thus for the annulus 



Zab{S) =Tre 



-5H„ 



Trg L -W2\ 



with q = e nS . As before, this can be decomposed into characters 

Z ab{5) =J2 U abXA(q), 

A 



(29) 
(30) 



but note that now the expression is linear. The non-negative integers n^ b give the operator 
content with the boundary conditions (aft): the lowest value of A with n^ b > gives the 
conformal weight of the bcc operator, and the others give conformal weights of the other 
allowed primary fields which may also sit at this point. 

On the other hand, the annulus partition function may be viewed, up to an overall rescal- 
ing, as the path integral for a CFT on a circle of unit circumference, being propagated 
for (imaginary) time <5 -1 . From this point of view, the partition function is no longer a 
trace, but rather the matrix element of e~ H ' 5 between boundary states: 



Z ab (5) = (a\e- H / s \b) 



(31) 



Note that H is the same hamiltonian that appears on the cylinder, and the boundary 
states lie in the Hilbert space of states on the circle. They can be decomposed into linear 
combinations of states in the representation spaces of the two Virasoro algebras, labelled 
by their lowest weights (A, A). 

How are these boundary states to be characterized? Recalling that on the cylinder L n oc 
f e mu T(u)du, and L n oc / e~ mu T(u)du, the conformal boundary condition implies that 
any boundary state \B) lies in the subspace satisfying 



L n \B) = L- n \B) 



(32) 
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This condition can be applied in each subspace. Taking n = in (132]) constrains A = A. 
It can then be shown that the solution of (132]) is unique within each subspace and has the 
following form. The subspace at level N of has dimension d&(N). Denote an ortho normal 
basis by \A,N;j), with 1 < j < d A (N), and the same basis for the representation space 
of Vir by |A, N;j). The solution to ( 1321) in this subspace is then 

oo d A (N) 



|A»=£ £ |A,iV;i)®|A,iV;i). (33) 

N=0 j=l 

These are called Ishibashi states. One way to understand this is to note that (132]) implies 
that 

(B\L n \B) = (B\i_ n \B) = (B\hatL n \B), 

where we have used L_ n = L n and assumed that the matrix elements are all real. This 

means that acting with the raising operators L n on \B) has exactly the same effect as the 
L n , so, starting with N = we must build up exactly the same state in the two spaces. 

Matrix elements of the translation operator along the cylinder between Ishibashi states 
are simple: 

((A'\e' &/S \A)) (34) 

oo d A ,(N') oo d A (N) 

= E E E E (&',N';f\(g)(A',N';j'\e~^^ Lo+Lo - c W 

N'=0 j'=l N=0 j=l 



\A,N;j)®\A,N;j) (35) 

oo d A (N) 

= 5 A , A Y, E e-^( A+N -^ = 5 A >AXA(e- A ^ s ) . (36) 

N=0 j=l 

Note that the characters which appear are related to those in (|30|) by the modular trans- 
formation 5*. 

The physical boundary states satisfying (13"U1) are linear combinations of these Ishibashi 
states: 

\a) = E((A|a) |A» . (37) 

A 

Equating the two different expressions (I30II31I) for Z a b, and using the modular trans- 
formation law for the characters and their linear independence, gives the (equivalent) 
conditions: 



n 



ab 



£<Sft<a|A'»«A'|&); (3? 



<a|A'»«A'|6> = Y, s £*&- ( 39 ) 

A 

The requirements that the right hand side of ([3*8]) should give a non-negative integer, 
and that the right hand side of (139]) should factorize in a and b, give highly nontrivial 
constraints on the allowed boundary states and their operator content. 



CFT and Statistical Mechanics 33 



, A 



// 



A' 

Figure 13: Argument illustrating the fusion rules. 

For the diagonal CFTs considered here (and for the nondiagonal minimal models) a 
complete solution is possible. Since the elements S^ of S are all non-negative, one may 
choose ((A|0) = (S^) . This defines a boundary state 

|6>=£(S A ) 1/2 |A», (40) 

A 

and a corresponding boundary condition such that n^ = <5ao- Then, for each A' ^ 0, one 
may define a boundary state 

«A|A')^,/(<?o A ) 1/2 - (41) 

From (|38|) . this gives n A , = <5a'a- For each allowed A' in the torus partition function, 
there is therefore a boundary state |A') satisfying (I38|l39p . However, there is a further 
requirement: 



n: 



cA q£ Qi 
A V^ J i °A'°A" 

A'A" 



J2^§^ (42) 



s, 



should be a non-negative integer. Remarkably, this combination of elements of S occurs 
in the Verlinde formula, which follows from considering consistency of the CFT on the 
torus. This states that the right hand side of (j4"2|) is equal to the fusion rule coefficient 
N A , A „. Since these are non- negative integers, the above ansatz for the boundary states is 
consistent. The appearance of the fusion rules in this context can be understood by the 
following argument, illustrated in Fig. [T3J Consider a very long strip. At 'time' t — > — oo 
the boundary conditions on both sides are those corresponding to 0, so that only states 
in the representation propagate. At time t\ we insert the bcc operator (0|A ; ) on one 
edge: the states A' then propagate. This can be thought of as the fusion of with A'. 
At some much later time we insert the bcc operator (0|A") on the other edge: by the 
same argument this should correspond to the fusion of A' and A", which gives all states 
A with A^, A „ = 1. But by definition, these are those with n A , A „ = 1. 

We conclude that, at least for the diagonal models, there is a bijection between the 
allowed primary fields in the bulk CFT and the allowed conformally invariant boundary 
conditions. For the minimal models, with a finite number of such primary fields, this 
correspondence has been followed through explicitly. 
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7.2.1 Example 

The simplest example is the diagonal c—\ unitary CFT corresponding to p — 4,p' = 3. 
The allowed values of the conformal weights are h — 0, ^A, and 
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from which one finds the allowed boundary states 



10} 
ll> 



"y3 1 )) + 72I2)) ~~ W^lie)) ! 
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The nonzero fusion rule coefficients of this CFT are 



iV o,o 



u, 16 



N k 



N'{ 

16 : 



JVl 



TVVi 

2'2 



A^i' 



16 

_ 1 

16 '2 



(43) 



(44) 
(45) 
(46) 



The c = I CFT is known to describe the continuum limit of the critical Ising model, in 

which spins s = ±1 are localized on the sites of a regular lattice. The above boundary 

conditions may be interpreted as the continuum limit of the lattice boundary conditions 

s — 1, s — — 1 and free (/), respectively. Note there is a symmetry of the fusion rules which 

means that one could equally well have reversed the first two. This shows, for example, 

that the for (//) boundary conditions the states with lowest weights (corresponding 

to the identity operator) and | (corresponding to the the magnetisation operator at the 

boundary) can propagate. Similarly, the scaling dimension of the (/| ± 1) bcc operator is 
j_ 

16' 

7.3 Boundary operators and SLE 



Let us now apply the above ideas to the A m models. There should be a set of conformal 
boundary states corresponding to the entries of first row (r, 1) of the Kac table, with 
1 < r < m. It is an educated guess (confirmed by exact calculations) that these in fact 
correspond to lattice boundary conditions where the heights on the boundary are fixed to 
be at a particular node r of the A m graph. What about the boundary condition changing 
operators? These are given by the fusion rules. In particular, since (suppressing the index 
s = l) 



K 



r.2 



Su 



>\r-r'\,l 1 

we see that the bcc operator between r and r ± 1, corresponding to a single cluster 
boundary intersecting the boundary of the domain, must be a (2, 1) operator in the Kac 
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tableo This makes complete sense: if we want to go from r\ to r 2 we must bring together 
at least \r± — r 2 \ cluster boundaries, showing that the leading bcc operator in this case is at 
(|t"i — T2I , 1), consistent once again with the fusion rules. If the bcc operators corresponding 
to a single curve are (2, 1) this means that the corresponding states satisfy 

(2L- 2 -(2/g)L 2 _ 1 )\<h,i) = 0. (47) 

We are now going to argue that ( 1471) is equivalent to the statement that the cluster 
boundary starting at this boundary point is described by SLE. In order to avoid being 
too abstract initially, we'll first show how the calculations of a particular observable agree 
in the two different formalisms. 

Let ( be a point in the upper half plane and let P(C) be the probability that the curve, 
starting at the origin, passes to the left of this point (of course it is not holomorphic). 
First we'll give the physicist's version of the SLE argument (assuming familiarity with 
Werner's lectures). We imagine making the exploration process for a small Loewner time 
5t, then continuing the process to infinity. Under the conformal mapping f$t{z) which 
removes the first part of the curve, we get a new curve with the same measure as the 
original one, but the point ( is mapped to fst(C)- But this will lie to the right of the 
new curve if and only if the original point lay to the right of the original curve. Also, by 
integrating the Loewner equation starting from fo(z) = z, we have approximately 

9 J\+ 

fst(z) « Z+ + y/K8B t , 

z 



at least for z ^$> 5t. Thus we can write down an equation 2 ' '• 

P(Q = E 



P U + 2 ~f + V^SB t ) 



SB t 



where E[. . .]sB t means an average over all realisations of the Brownian motion up to time 
St. Expanding the right hand side to 0(St), and remembering that E[5S t ] = and 
E[(8B t ) 2 ] = 8t, we find (with ( = x + iy), the linear PDE 

( 2x d 2y d k d 2 \ „. , ..„. 

' ] P(x,y) = 0. (48) 



\x 2 + y 2 dx x 2 + y 2 dy 2 dx 2 

By scale invariance, P(x,y) depends in fact only on the ratio y/x, and therefore this can 
be reduced to a second order ODE, whose solution, with appropriate boundary condi- 
tions, can be expressed in terms of hypergeometric functions (and is known as Schramm's 
formula. ) 



24 If instead of the dilute lattice model we consider the dense phase, which corresponds, eg to the 
boundaries of the FK clusters in the Potts model, then r and s get interchanged for a given central 
charge c, and the bcc operator then lies at (1,2). 

25 Some physicists will recognise this as the reverse Fokker-Planck equation. 
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Now let us give the CFT derivation. In terms of correlation functions, P can be expressed 
as 

p _ (02,l(O)$(C, 002,1 (°0)) 

(02,l(O)02,l(oo)) 

The denominator is just the partition function restricted to there being a cluster boundary 
from to infinity. $ is an 'indicator operator' which takes the values or 1 depending 
on whether the curve passes to the right (respectively left) of (. Since P is a probability 
it is dimensionless, so $ has zero conformal dimensions and transforms trivially. 

Now suppose we insert J c (2T(z)/z)(dz/2iri) + cc. into the correlation function in the 
numerator, where C is a semicircular contour surrounding the origin but not (. Using the 
OPE of T with 02,i, this gives 

2L_ 2 2il (O) = (2/^)^0 li2 (x)| x=o . 

Using translation invariance, this derivative can be made to act equivalently on the x- 
coordinate of (. On the other hand, we can also distort the contour to wrap around £, 
where it simply shifts the argument of $. The result is that we get exactly the same PDE 
as in (|48p . with the identification 

g = A/K. 

Of course this was just one example. Let us see how to proceed more generally. In 
radial quantisation, the insertion of the bcc field 02,i(O) gives a state \<f>2,i)- Under the 
infinitesimal mapping f dt we get the state 

(l - (2L_ 2 dt + L- X <s/tidBt)) |0 2 ,i) , 

or, over a finite time, a time-ordered exponential 

T exp (- J (L_ 2 di' + L- iy /KdBt,)) \fa,i) . (49) 

The conformal invariance property of the measure on the curve then implies that, when 
averaged over dB t >, this is in fact independent of t. Expanding to 0(t) we then again 
find ( H71) with g = 4/k. Since this is a property of the state, it implies an equivalence 
between the two approaches for all correlation functions involving 2) i(O), not just the one 
considered earlier. Moreover if we replace y/ndB t by some more general random driving 
function dW t , and expand (T49|) to any finite order in t using the Virasoro algebra and 
the null state condition, we can determine all moments of W t and conclude that it must 
indeed be rescaled Brownian motion. 

Of course the steps we used to arrive at this result in CFT are far less rigorous than the 
methods of SLE. However, CFT is more powerful in the sense that many other similar 
result can be conjectured which, at present, seem to be beyond the techniques of SLE. 
This is part of an ongoing symbiosis between the disciplines of theoretical physics and 
mathematics which, one hopes, will continue. 
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The basic reference for CFT is the 'big yellow book': Conformal Field Theory by P. di 
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Statistical Field Theory by C. Itzykson and J.-M. Drouffe (Cambridge University Press, 
1989.) A gentler introduction is provided in the 1988 les Houches lectures by P. Ginsparg 
and J. Cardy in Fields, Strings and Critical Phenomena, E. Brezin and J. Zinn- Justin, 
eds. (North-Holland, 1990.) Other specific pedagogical references are given below. 
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Press, 1996.) 
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6.1 V. Pasquier, Nucl. Phys. B 285, 1 62, 1986, and J. Phys. A 20, L1229, 1987; I. Ros- 
tov, Nucl. Phys. B 376, 539, 1992; |arXiv:hep-th/9112059| The version in these lec- 
tures is discussed in J. Cardy, J.Phys. A 40, 1427 (2007); [arXrV:math-ph/0610030[ 

6.2 The basic reference is B. Nienhuis, in Phase transitions and critical phenomena, 
vol. 11, p. 1, C. Domb and J.L. Lebowitz, eds. (Academic, 1987.) A slightly 
different approach is explained in J. Kondev, Phys. Rev. Lett. 78, 4320, 1997. 

7.1,7.2 For reviews see V.B. Petkova and J.-B. Zuber, Conformal boundary conditions and 
what they teach us, lectures given at the Summer School and Conference on Non- 
perturbative Quantum Field Theoretic Methods and their Applications, August 



2000, Budapest, Hungary, arXiv:hep-th/0103007, and J. Cardy, Boundary Con- 



formal Field Theory in Encyclopedia of Mathematical Physics, (Elsevier, 2006); 



arXiv:hep-th/0411189. 



7.3 For reviews on the connection between SLE and CFT, see J. Cardy, SLE for The- 
oretical Physicists, Ann. Phys. 318, 81-118 (2005); |arXiv:cond-mat/0503313~ 



and M. Bauer and D. Bernard, 2D growth processes: SLE and Loewner chains, 
Phys.Rept. 432, 115-221 (2006); |arXiv:math-p h/0602049. 



